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Abstract 

The  spatial  scale  on  which  microbial  communities  respond  to  plant  invasions  may 
provide  important  clues  as  to  the  nature  of  potential  invader-microbe  interactions. 
Lespecleza  cuneata  (Dum.  Cours.)  G.  Don  is  an  invasive  legume  that  may  benefit  from 
associations  with  mycorrhizal  fungi;  however,  it  has  also  been  suggested  that  the  plant  is 
allelopathic  and  may  alter  the  soil  chemistry  of  invaded  sites  through  secondary 
metabolites  in  its  root  exudates  or  litter.  Thus,  L.  cuneata  invasion  may  interact  with  soil 
microorganisms  on  a  variety  of  scales.  We  investigated  L.  cuneata- related  changes  to  soil 
bacterial  and  fungal  communities  at  two  spatial  scales  using  multiple  sites  from  across  its 
invaded  N.  American  range.  Using  whole  community  DNA  fingerprinting,  we 
characterized  microbial  community  variation  at  the  scale  of  entire  invaded  sites  and  at  the 
scale  of  individual  plants.  Based  on  permutational  multivariate  analysis  of  variance,  soil 
bacterial  communities  in  heavily  invaded  sites  were  significantly  different  from  those  of 
uninvaded  sites,  but  bacteria  did  not  show  any  evidence  of  responding  at  very  local  scales 
around  individual  plants.  In  contrast,  soil  fungi  did  not  change  significantly  at  the  scale  of 
entire  sites,  but  there  were  significant  differences  between  fungal  communities  of  native 
versus  exotic  plants  within  particular  sites.  The  differential  scaling  of  bacterial  and  fungal 
responses  indicates  that  L.  cuneata  interacts  differently  with  soil  bacteria  and  soil  fungi, 
and  these  microorganisms  may  play  very  different  roles  in  the  invasion  process  of  this 
plant. 
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1.  Introduction 


A  growing  body  of  research  indicates  that  soil  microorganisms  can  both  respond 
to  and  strongly  impact  plant  communities,  including  those  affected  by  invasive  plants 
(Klironomos  2002,  Kourtev,  et  al.  2003,  Reinhart,  et  al.  2003,  Callaway,  et  al.  2004, 
Reinhart  and  Callaway  2004,  van  der  Putten,  et  al.  2007b,  Engelkes,  et  al.  2008); 
reviewed  by  (Wardle,  et  al.  2004,  Wolfe  and  Klironomos  2005,  Reinhart  and  Callaway 
2006,  van  der  Putten,  et  al.  2007a).  Plant-microbe  interactions  can  influence  the 
competitive  dynamic  of  plant  communities  (Bever  1994,  Bever  2003),  explain  variation 
in  the  local  abundance  of  plant  species  (Klironomos  2002),  and  may  facilitate  range 
expansion  and  invasion  (Engelkes,  et  al.  2008).  Our  understanding  of  the  specific 
microbial  mechanisms  at  play  in  plant  community  dynamics  is  limited  because  much 
previous  work  has  considered  microbial  communities  from  a  purely  functional,  "black 
box"  perspective;  that  is,  a  particular  microbial  species  pool  either  net  positively  or  net 
negatively  influences  plant  growth  (examples:  (Reinhart,  et  al.  2003,  Agrawal,  et  al. 

2005,  Zou,  et  al.  2006)).  Furthermore,  attempts  to  correlate  plant  invasion  to  specific 
changes  in  microbial  community  composition  have  met  with  mixed  results  (Kourtev,  et 
al.  2003,  Callaway,  et  al.  2004,  Kardol,  et  al.  2007,  Lorenzo,  et  al.  2010),  and  much  is 
still  unknown,  in  particular  about  the  specific  microbial  communities  that  respond  to 
plant  invasions. 

The  question  of  scale  is  one  potential  hurdle  to  a  deeper  understanding  of 
microbial  ecology's  role  in  plant  invasions,  because  different  interactions  can  operate  on  a 
number  of  scales.  Microbes  can  interact  directly  with  invaders  on  very  fine  scales  through 
a  number  of  symbiotic  interactions  with  roots:  e.g.  mycorrhizas,  mutualistic  bacteria,  or 
known  pathogens  (Klironomos  2002,  Klironomos  2003,  Fumanal,  et  al.  2006,  Reinhart 
and  Callaway  2006,  van  der  Putten,  et  al.  2007b,  Pringle,  et  al.  2009).  It  is  well  known 
that  the  rhizosphere  is  a  “hot  spot”  of  interaction  that  can  influence  microbial  activity  and 
plant  ecological  performance  (Kent  and  Triplett  2002,  Berg  and  Smalla  2009).  However, 
the  soil  milieu  is  broader  than  the  rhizosphere,  and  understanding  the  microbial  ecology 
of  some  plant  invasions  may  require  observations  at  larger  scales,  where  interactions  may 
be  more  indirect. 

Several  recent  investigators  have  proposed  that  invasive  plant  impacts  may  be 
more  diffuse,  affecting  microbial  communities  that  are  not  part  of  the  rhizosphere  (Wolfe 
and  Klironomos  2005,  van  der  Putten,  et  al.  2007a).  For  instance,  through  inputs  of  litter 
to  the  soil,  invasive  plants  may  alter  the  decomposition  rate  and  nutrient  cycling  rates  of 
the  system  (Ehrenfeld  2003,  Heneghan,  et  al.  2006,  Zou,  et  al.  2006),  with  consequences 
for  nutrient  availability  and  competitive  fitness  of  native  competitors.  Many  plant 
invaders  have  been  found  to  alter  the  soil  characteristics,  such  as  fungakbacterial  biomass 
ratios  (Kourtev,  et  al.  2002,  Kourtev,  et  al.  2003),  and  total  nitrogen,  carbon  or 
phosphorus  pools  (Ehrenfeld  2003,  Chapuis-Lardy,  et  al.  2006,  Heneghan,  et  al.  2006, 
Lorenzo,  et  al.  2010),  and  these  observations  suggest  that  invasive  plants  may  influence 
soil  ecology  on  scales  larger  than  that  of  their  own  rhizosphere  and  over  longer  time 
periods  than  the  lifetime  of  individual  plants.  The  “novel  weapons”  hypothesis  suggests 
that  invasive  plants  may  also  introduce  new  chemical  compounds  into  the  soil 


(Cappuccino  and  Amason  2006,  Callaway,  et  al.  2008),  and  the  effects  of  these 
compounds  will  not  necessarily  be  restricted  to  rhizosphere  communities.  For  example, 
there  is  evidence  that  allelochemicals  secreted  by  invasive  garlic  mustard  ( Alliaria 
petiolata )  (Stinson,  et  al.  2006,  Wolfe,  et  al.  2008,  Lankau  2010)  and  spotted  knapweed 
(Cenaurea  maculosa)  (Callaway,  et  al.  2004,  Broz,  et  al.  2007)  can  reduce  the  abundance 
and  diversity  of  soil  fungi  that  are  important  mutualists  of  native  competitor  plants  (the 
“degraded  mutualisms”  hypothesis  (Vogelsang  and  Bever  2009)).  Invader  impacts 
extending  beyond  the  rhizosphere  can  have  important  implications  for  the  invaded 
habitats.  Thus,  while  rhizosphere  interactions  are  very  important,  it  is  also  necessary  to 
understand  how  invasive  plants  can  impact  the  soil  at  larger  spatial  scales  where  we  can 
find  the  signature  of  other  important  interactions. 

The  Chinese  bush  clover  ( Lespedeza  cuneata  (Dum.  Cours.)  G.  Don),  often  called 
Sericea  Lespedeza,  is  one  of  about  45  Lespedeza  species  found  across  Asia,  Australia, 
and  North  America.  These  legumes  are  tolerant  of  degraded  and  nutrient-poor  soils,  and 
they  have  been  shown  to  associate  with  a  broad  diversity  of  rhizobacteria  (Gu,  et  al. 
2007).  L.  cuneata  is  a  long-lived  perennial  introduced  from  Japan  to  the  United  States  in 
the  late  1800s  (Guernsey  1970),  and  it  is  now  well  established  as  an  invasive  weed  of  old 
fields,  prairies  and  forest  openings  across  the  southern  U.S.  (Miller  2003).  L.  cuneata  has 
a  high  seed  production  rate  and  a  high  population  growth  rate  compared  with  its  North 
American  congeners  (Schutzenhofer,  et  al.  2009,  Woods,  et  al.  2009).  It  can  form  dense 
stands  in  its  invaded  range,  shading  out  native  grasses  and  herbs  (Brandon,  et  al.  2004). 
There  is  some  evidence  that  its  residues  possess  phytotoxic  or  allelopathic  chemicals 
(Langdale  and  Giddens  1967,  Kalburtji,  et  al.  2001),  and  its  residues  and  root  exudates 
have  been  found  to  inhibit  seed  germination  and/or  growth  of  various  grass  species 
(Kalburtji  and  Mosjidis  1992,  Kalburtji  and  Mosjidis  1993b,  Kalburtji  and  Mosjidis 
1993a).  In  addition  to  its  rhizobacterial  symbionts,  it  has  been  found  to  associate  with 
various  mycorrhizal  fungi  in  what  has  been  termed  a  nitrogen-fixing,  tri-partite  symbiosis 
(Lynd  and  Ansman  1993). 

Due  to  its  symbiotic  associations,  L.  cuneata  invasions  may  interact  with  soil 
microbial  communities  of  invaded  systems  at  very  local  scales  (i.e.  that  of  the  individual 
plant).  In  addition,  the  potential  inhibitory  effects  of  L.  cuenata  root  exudates  and 
residues  suggest  that  this  invader  may  affect  the  soil  system  at  broader  scales,  particularly 
where  it  achieves  high  population  densities.  In  this  study,  we  investigate  variation  in  soil 
bacterial  and  fungal  community  composition  associated  with  L.  cuneata  invasions  at  two 
spatial  scales:  1)  in  localized  soil  zones  around  individual  plants  (hereafter  referred  to  as 
plant  “neighborhoods”  and  2)  that  of  whole  sites,  including  sites  with  sparse  and  those 
with  dense  populations  of  L.  cuneata.  Because  its  broad  geographic  distribution  in  the 
U.S.  means  that  L.  cuneata  may  interact  with  many  different  combinations  of  local  soil 
microbiota,  we  seek  patterns  that  are  consistent  across  the  invaded  range,  and  our 
sampling  across  multiple  widely  separated  locations  represents  a  third,  geographic  scale 
of  microbial  community  variation. 


2.  Materials  and  Methods 


2.1  Study  design.  To  evaluate  potential  impacts  of  Lespedeza  cuneata  invasion  intensity 
on  soil  microbial  communities,  field  sampling  was  conducted  according  to  a  blocked 
sampling  design,  with  the  blocks  consisting  of  five  different  U.S.  military  installations 
(hereafter  referred  to  as  regions)  distributed  throughout  the  invaded  range  of  L.  cuneata 
in  the  southern  and  midwestern  USA  (Table  1).  Each  installation  was  sampled  once 
between  late  May  and  late  July  2009.  At  each  installation,  interviews  with  the  natural 
resource  managers  were  used  to  identify  three  sampling  sites  (Fig.  1)  differing  in  invasion 
intensity:  1)  “uninvaded”  sites  contained  no  L.  cuneata  individuals;  2)  “lightly  invaded” 
sites  contained  L.  cuneata  as  a  low  abundance  member  of  the  herbaceous  community; 
and  3)  “heavily  invaded”  sites  were  dominated  by  L.  cuneata.  These  invasion  levels 
constitute  an  experimental  treatment  at  the  whole-site  scale,  replicated  across  all  five 
installations  (regions).  Sites  ranged  in  area  from  ~40  m“  to  -4500  nr,  and  the  distance 
between  sites  ranged  from  100-200m  (Ft.  Riley),  to  l-3km  (Fts.  Benning,  Campbell,  and 
Feonard  Wood),  to  10-25km  (Ft.  Stewart).  Characteristics  of  sites,  such  as  soil  type, 
habitat,  and  the  identity  of  other  leguminous  plant  species,  differed  broadly  across  the 
installations  (Table  1). 


Table  1.  Site  characteristics  for  five  military  installations  sampled  in  2009. 


Region 

Location 

INVASION 

LEVEL 

Habitat 

Soil* 

Exotic  plants 

Native  plants 

uninvaded 

forest 

Fuquay  loamy  sandb 

L.  stuevei 

L.  violacea 

Stewart 

E  Georgia 

light 

open 

field 

Fuquay  loamy  sandb 

L.  cuneata 

C.  fasciculata 

D.  rotundifolium 

heavy 

open 

field 

Mascotte  fine  sandc 

L.  cuneata 

M.  lupulina 

C.  fasciculata 

uninvaded 

forest 

Fuquay  loamy  sandb 

L.  hirta 

L.  repens 

L.  virginica 

Benning 

W  Georgia 

light 

forest 

Cowarts  fine  sandy  loamb 

L.  cuneata 

C.  fasciculata 

L.  virginica 

heavy 

open 

field 

loamy  Udorthents* 

L.  cuneata 

T.  campestre 

uninvaded 

open 

field 

Pickwick  silt  loamb 

L.  procumbens 

L.  virginica 

Campbell 

SW  Kentucky 

light 

forest 

Moundville  silt  loamb 

L.  cuneata 

C.  fasciculata 

L.  virginica 

heavy 

open 

field 

Dickson  silt  loamb 

L.  cuneata 

K.  striata 

C.  fasciculata 

uninvaded 

forest 

Union  silt  loamd 

D.  paniculatum 

L.  virginica 

Leonard 

Wood 

Central 

Missouri 

light 

forest 

Union  silt  loamd 

L.  cuneata 

C.  fasciculata 

L.  virginica 

heavy 

open 

field 

Sturkie  silt  loam6 

L.  cuneata 

C.  varia 

uninvaded 

prairie 

Smolan  silt  loam6 

L.  capitata 

Riley 

NE  Kansas 

light 

prairie 

Geary  silt  loam6 

L.  cuneata 

C.  fasciculata 

L.  capitata 

heavy 

open 

field 

Geary  silt  loam6 

L.  cuneata 

C.  fasciculata 

L.  capitata 

a  obtained  from  National  Cooperative  Soil  Survey  (NCSC  Data,  http://ssldata.nrcs.usda.com) 

b  f  Soil  orders: b  Ultisol, c  Spodosol, d  Alfisol, e  Mollisol, f  Udorthents  =  highly  disturbed,  with  little  similarity  to  the  original  soil 


At  each  site  in  each  installation,  four  individual  plants  of  at  least  two  different 
leguminous  species  were  collected  for  characterization  of  microbial  communities  in  the 
soil  of  their  root  balls.  These  samples  were  used  to  determine  if  different  plants  alter  the 


soil  microbial  community  composition  of  their  local  soil  “neighborhoods.”  L.  cuneata 
was  sampled  at  all  lightly  and  heavily  invaded  sites,  along  with  at  least  one  native  legume 
species.  Only  multi-stemmed  L.  cuneata  were  selected  for  sampling,  in  order  to  survey 
the  older  individuals  at  each  site.  An  additional  exotic  legume  was  sampled  at  most  of  the 
heavily  invaded  sites  (Table  1),  and  all  other  plants  sampled  were  native  species.  Root 
ball  soils  were  collected  by  digging  -10-15  cm  around  the  stem  of  a  plant,  to  a  depth  of 
-15  cm,  and  then  extracting  the  entire  intact  root  ball  into  a  plastic  bag,  which  was  placed 
on  ice  for  transport  to  the  laboratory.  At  the  lab,  plants  were  vigorously  shaken  to  release 
the  soil  from  the  root  ball,  and  this  soil  was  sub-sampled  into  a  sterile  15-ml  centrifuge 
tube  and  frozen  until  further  processing  (typically  less  than  two  weeks).  At  all  lightly  and 
heavily  invaded  sites,  bulk  soil  was  collected  from  a  plant-free  location  by  digging  a  hole 
with  the  previously  described  dimensions.  Bulk  soil  was  otherwise  treated  as  described 
for  root  ball  soil  samples  in  all  respects.  All  samples  were  geo-referenced  using  a  hand¬ 
held  GPS,  and  these  coordinates  were  used  to  obtain  soil  data  (order  and  series)  for  each 
site  using  the  National  Resources  Conservation  Service’s  Web  Soil  Survey  (Soil  Survey 
Staff). 


The  study  design  incorporates  three  spatial  scales  of  investigation.  The  five 
regions  served  as  blocks  (i.e.  replicates)  of  the  three-level  treatment,  invasion  level 
(“uninvaded”,  “lightly  invaded”,  and  “heavily  invaded”),  for  a  total  of  15  sites  distributed 
on  a  large  geographical  scale  throughout  the  invaded  range  of  L.  cuneata.  Site-to-site 
comparisons  (e.g  invasion  level)  characterize  changes  on  a  smaller  (regional)  scale.  For 
analyses  at  these  scales,  samples  collected  from  different  plants  or  bulk  soil  were  treated 
as  subsamples  (see  Data  Analysis),  because  they  did  not  constitute  independent 
observations  of  microbial  communities  at  any  particular  site.  However,  the  design  allows 
for  valid  comparisons  of  small-scale  differences  between  different  plants  at  the  same  site 
using  the  (n  =  4)  replicates  of  each  plant  species  sampled.  These  samples  comprised  a  16- 
level  treatment  (plant)  that  was  nested  within  the  15  main  sites.  Many  plants  were 
restricted  to  a  single  facility,  and  some  plants  were  only  sampled  at  one  invasion  level  at 
a  given  facility  (Table  1).  All  plants  sampled  were  legumes,  and  all  but  five  species  were 
natives.  The  different  plant  species  sampled  (Table  1)  included:  Lespedeza  cuneata,  L. 
captitata ,  L.  hirta ,  L.  procumbens,  L.  repens ,  L.  stuevei,  L.  violacea,  L.  virginica, 
Chamaecrista  fasciculata,  Coronilla  varia,  Desmodium  paniculatum ,  D.  rotundifolium, 
Kummerowia  striata,  Meclicago  lupulina,  and  Trifolium  campestre.  For  data  analysis 
purposes,  bulk  soil  cores  were  considered  to  comprise  a  single  “species”  in  the  (nested) 
treatment  plant. 

2.2  DNA  extraction  and  ARISA.  Soil  samples  were  lyophilized  for  24-48  h  until 
completely  dry  and  then  homogenized  by  hand  with  a  sterile  metal  pestle.  DNA  was 
extracted  from  300  mg  of  lyophilized  soil  from  each  sample  using  the  FastDNA®  SPIN 
Kit  for  Soil  (MP  Biomedicals,  LLC,  Solon,  OH,  USA)  according  to  the  manufacturer’s 
instructions.  The  DNA  was  further  purified  to  remove  potential  PCR-inhibiting 
contaminants  through  a  15 -minute  incubation  at  65°C  with  CTAB  (cetyl  trimethyl 
ammonium  bromide;  final  concentration:  1%  CTAB,  0.7M  NaCl),  followed  by  extraction 
with  24: 1  chloroforimisoamyl  alcohol  and  precipitation  with  EtOH.  Pellets  were  washed 
twice  with  70%  EtOH  and  dissolved  in  IX  TE  buffer. 


Microbial  community  composition  was  assessed  using  ARISA,  or  Automated 
Ribosomal  Intergenic  Spacer  Analysis  (Fisher  and  Triplett  1999).  ARISA  profiles  were 
generated  for  both  bacterial  and  fungal  communities.  Bacterial  ARISA  targeted  the  16S- 
ITS-23S  regions  of  bacterial  rm  operons,  using  the  primer  set  of  Borneman  and  Triplett 
(Borneman  and  Triplett  1997):  1406f  (5’-TGyACACACCGCCCGT-3’)  and  23Sr  (5’- 
GGGTTbCCCCATTCrG-3’).  Fungal  ARISA  targeted  the  18S-ITS1-5.8S-ITS2-28S 
regions  of  fungal  nuclear  ribosomal  DNA  using  the  primer  set  of  Ranjard  and  colleagues 
(Ranjard,  et  al.  2001):  2234C  (5 ’ - GTTTC C GT AGGTG A ACCT GC- 3 ’ )  and  3126T  (5’- 
ATATGCTTAAGTTCAGCGGGT-3’).  Fluorochromes  were  added  to  the  5’  ends  of  the 
primers  1406f  (6-FAM)  and  3126T  (HEX)  so  that  PCR  products  could  be  visualized 
during  capillary  gel  electrophoresis.  All  PCR  reactions  contained  IX  buffer  (50  mM  Tris, 
pH  8.3),  250  mg/ml  bovine  serum  albumin,  3.0  mM  MgCl2,  250  mM  of  each  dNTP,  400 
nM  of  the  appropriate  forward  and  reverse  primers,  1.25  U  of  Taq  polymerase  (Promega, 
Madison,  WI,  USA),  and  20  ng  of  template  DNA  in  a  final  volume  of  25  ml.  PCR 
conditions  consisted  of  a  2-min.  initial  denaturation  at  94°C,  followed  by  30  cycles  of 
94°C  for  35s,  55°C  for  45s,  72°C  for  2  min.,  and  a  final  extension  at  72°C  for  2  min. 
Capillary  gel  electrophoresis  to  resolve  ARISA  fragments  was  performed  on  an  ABI 
Prism  3730x1  Analyzer  by  the  W.M  Keck  Center  for  Comparative  and  Functional 
Genomics,  University  of  Illinois  at  Urbana-Champaign. 

ARISA  profiles  were  analyzed  using  the  software  GeneMarker  v.  1.85 
(SoftGenetics,  LLC,  State  College,  PA,  USA),  which  converted  the  fluorescence  data  into 
electropherograms  and  estimated  the  sizes  (in  bps)  of  ARISA  fragments  based  on  an 
internal  size  standard  (MM-1000-ROX,  BioVentures,  Inc.,  Murfreesboro,  TN,  USA). 
ARISA  fragments  of  similar  sizes  (within  ~  1-2  bp  of  each  other)  were  considered  to 
represent  the  same  microbial  taxa  across  all  samples,  and  they  were  therefore  “binned” 
together  to  create  a  set  of  Operational  Taxonomic  Units  (OTUs)  for  this  study.  ARISA 
fragment  bins  were  generated  automatically  with  the  GeneMarker  software,  and  the 
resulting  set  of  bins  was  manually  corrected  to  remove  any  overlap  between  bins,  delete 
bins  generated  by  spurious  peaks,  and  create  additional  bins  to  accommodate  small  peaks 
missed  by  the  software.  Only  peaks  with  signal  intensities  >100  relative  fluorescence 
units  and  belonging  to  the  following  size  ranges  were  considered  valid:  400-1000bp  for 
bacteria  (Fisher  and  Triplett  1999),  and  300-1000bp  for  fungi  (Ranjard,  et  al.  2001). 
ARISA  fragment  intensity  was  determined  from  the  area  of  each  peak,  and  the  intensity 
value  and  bin  indentifications  were  exported  for  all  peaks  from  all  samples  to  create  a 
sample -by-OTU  data  table. 

2.3  Data  analysis.  Bacterial  and  fungal  communities  were  analyzed  separately  using 
multivariate  techniques.  ARISA  data  were  first  Hellinger-transformed,  which  consists  of: 
1)  dividing  each  ARISA  peak  area  by  the  sum  of  all  peak  areas  in  the  sample,  and  2) 
taking  the  square  root  of  each  value  obtained  in  step  1  (Legendre  and  Gallagher  2001). 
The  Hellinger  transformation  converts  the  raw  intensity  of  each  peak  to  a  relative 
intensity,  which  controls  for  any  run-to-run  variability  in  signal  strength  that  may  arise 
during  capillary  electrophoresis.  The  Hellinger-transformed  data  were  then  converted  into 


two  distance  matrices  (i.e.  one  for  bacteria,  one  for  fungi)  using  the  Bray-Curtis  index  of 
community  similarity  (Legendre  and  Legendre  1998). 

Permutational  Multivariate  Analysis  of  Variance  (PerMANOVA)  (McArdle  and 
Anderson  2001)  was  used  to  test  hypotheses  regarding  regions  and  invasion  level 
effects.  This  procedure  partitions  the  sums  of  squares  of  multivariate  data  sets  in  a 
manner  similar  to  that  of  a  standard  Analysis  of  Variance,  except  that  it  can  use  any 
ecologically-relevant  distance  metric  between  sample  points  to  calculate  the  sums  of 
squares  (McArdle  and  Anderson  2001).  The  resulting  “pseudo-F”  statistic  is  then 
compared  to  the  (upper)  tail  of  a  null  distribution  of  pseudo-F’ s  generated  by  random 
permutation  (shuffling  the  raw  data  between  treatments).  If  the  original  pseudo-F  is 
greater  than  a  given  proportion  (p  =  a)  of  the  pseudo-F’ s  in  the  null  distribution,  then  the 
null  hypothesis  can  be  rejected  with  a  type  I  error  rate  of  a.  We  used  an  implementation 
of  PerMANOVA  that  involved  sequential  removal  of  sums  of  squares  attributable  to  each 
term  in  the  model,  and  thus  the  order  in  which  terms  are  specified  in  the  model  influences 
the  significance  level  of  the  terms.  This  allowed  us  to  specify  models  controlling  for  any 
potential  "nuisance"  variation  due  to  site-level  differences  in  soil  texture  and  habitat 
(Table  1).  For  site-scale  analyses,  we  chose  the  most  conservative  model  in  which 
invasion  level  was  entered  last,  meaning  that  we  were  testing  for  the  differences  due  to 
INVASION  LEVEL  after  removing  any  differences  due  to  REGIONS,  SOIL  TEXTURE,  and 
HABITAT. 

The  null  distribution  for  PerMANOVA  tests  for  site-scale  effects  was  generated 
using  a  restricted  permutation  scheme  that  maintained  the  grouping  of  subsamples 
collected  from  the  same  site  (i.e.  different  plant  and  bulk  soil  samples).  This  was 
necessary  because  these  samples  are  properly  considered  pseudoreplicates  at  the  scale  of 
whole  site,  and  treating  them  as  full  replicates  would  artificially  inflate  the  error  degrees 
of  freedom.  The  stratified  randomization  scheme  preserved  the  true  level  of  replication 
for  the  hypothesis  test  of  regions  and  invasion  level  (5  regions  x  3  invasion  levels 
=  15),  while  allowing  us  to  use  data  from  the  full  sample  set. 

To  test  the  hypothesis  that  different  species  of  plants  differentially  alter  the 
microbial  communities  of  their  own  soil  neighborhoods,  PerMANOVA  tests  were  carried 
out  separately  for  each  site.  We  did  not  use  a  stratified  randomization  scheme  for  these 
tests,  under  the  assumption  that  all  (n  =  4)  samples  from  each  plant  constituted 
independent  observations  of  soil  neighborhood  communities  for  the  various  sampled 
species  at  each  site.  Because  we  were  primarily  interested  in  plant  neighborhood  effects 
due  to  L.  cuneata,  these  analyses  were  only  conducted  for  the  10  sites  at  which  L. 
cuneata  was  present  (i.e.  lightly  and  heavily  invaded  sites;  Table  1). 

For  all  significant  overall  PerMANOVA  tests  (regions,  invasion  level,  plant, 
etc.),  post  hoc  pairwise  PerMANOVA  tests  were  performed  to  identify  significant 
difference  among  the  treatment  levels.  All  post  hoc  significance  tests  used  the  Bonferroni 
correction  to  preserve  the  overall  type  I  error  rate  of  a  =  0.05.  For  visualization  of  the 
data,  Nonmetric  multidimensional  scaling  (NMDS)  was  also  applied  to  all  samples,  but 
only  the  centroids  and  standard  deviations  of  the  different  treatment  groups  are  displayed 


on  the  plots.  NMDS  is  an  ordination  technique  that  plots  multivariate  data  such  that  the 
rank-order  of  distances  between  the  plot  points  (representing  the  samples)  are 
monotonically  related  to  the  pairwise  distances  between  samples  in  the  original  (ARISA) 
distance  matrix.  This  technique  has  been  previously  recommended  for  analysis  of 
microbial  community  data  exhibiting  high  complexity  and/or  high  beta  diversity 
(Culman,  et  al.  2008).  The  Bray-Curtis  distance  matrices  of  Hellinger-transformed 
ARISA  data  were  used  as  inputs  for  NMDS,  and  two-  or  three-dimensional  solutions 
were  sought  for  ease  of  visualization.  Initial  configurations  for  NMDS  were  derived  from 
Principal  Coordinates  Analysis  of  the  data,  and  all  NMDS  solutions  converged  on 
optimal  solutions  in  less  than  50  iterations.  Final  stress  was  calculated  as  Kruskal's  stress 
formula  1  (Venables  and  Ripley  2002). 

All  data  analyses  were  performed  in  the  R  statistical  environment  (R 
Development  Core  Team  2005)  using  functions  found  in  the  packages  "MASS" 

(Venables  and  Ripley  2002)  and  "vegan"  (Oksanen,  et  al.  2009).  Hellinger  transformation 
used  the  function  decostand  (method="hellinger" ) ,  PerMANOVA  used  the 
function  adonis  ( )  with  1000  permutations  of  the  data,  and  NMDS  used  the  function 
isoMDS  ()  . 

3.  Results 

Bacterial  communities  showed  significant  variation  due  to  different  invasion 
levels  of  L.  cuneata  even  after  the  effects  of  regions,  soil  texture,  and  habitat  had 
been  removed,  but  fungal  communities  did  not  (Table  2).  Pairwise  contrasts  revealed  that 
bacterial  communities  from  heavily  invaded  sites  were  significantly  different  from 
uninvaded  sites  (Fig.  2). 


Table  2.  Significance  levels  of  large-scale  effects*1 


Model  terms 

Bacteria 

p-value 

Fungi 

REGIONS 

0.001  *** 

0.001  *** 

SOIL  TEXTURE 

0.43 

0.02b  * 

HABITAT 

0.13 

0.08b 

INVASION  LEVEL 

0.01  * 

0.08 

“Significance  level  determined  by  permutational  MANOVA 
Switching  the  order  of  these  two  terms  in  the  model  for  fungi  also  reverses 
their  significance  levels 
*  0.01  <  p  <  0.05;  **  0.001  <  p  <  0.01 


Both  bacterial  and  fungal  community  composition  varied  significantly  across 
regions  (p  >  0.001  for  bacteria  and  fungi),  but  they  displayed  different  patterns  of 
geographic  variation  across  the  various  military  installations  (Fig.  3).  For  bacteria, 
pairwise  contrasts  indicated  that  the  communities  of  Ft.  Riley  (KS)  were  different  from 
those  of  Ft.  Leonard  Wood  (MO),  with  the  communities  of  Ft.  Campbell  (KY)  and  the 
two  GA  facilities  intermediate  to  these  two  end  groups  (Fig.  3A-B).  In  contrast,  fungal 
communities  at  Fts.  Stewart  (GA)  and  Riley  were  deemed  to  be  different  from  those  of 
Ft.  Campbell,  with  those  of  the  remaining  facilities  being  statistically  indistinguishable 
from  either  of  these  two  groups  (Fig.  3C-D). 


Neither  habitat  nor  soil  texture  was  related  to  significant  differences  in 
bacterial  communities  (Table  2).  For  fungi,  inclusion  of  habitat  and  SOIL  texture  in 
the  PerMANOVA  model  always  considered  the  first  term  added  to  be  significant,  but  not 
the  second  term  (e.g.  Table  2  displays  a  model  with  SOIL  texture  added  first).  Thus,  at 
least  one  of  these  factors  was  a  significant  cause  of  variation  in  soil  fungal  communities 
(Fig.  4).  Pairwise  contrasts  showed  that  fungi  from  forested  habitats  were  different  from 
those  of  prairie  and  open  field  habitats  (Fig.  4A).  Soil  texture  divided  fungal  communities 
into  three  principal  groups:  those  from  sandy  loam  and  loamy  sand  soils,  those  from 
sandy  soils,  and  those  from  loam  soils  (Fig.  4B).  Fungal  communities  on  silty  loam  soils 
were  not  significantly  different  than  any  of  these  groups,  although  they  were 
distinguishable  from  the  loamy  sand  communities  (Fig.  4B).  Fungal  communities  of 
heavily  invaded  sites  were  shifted  along  NMDS  axis  1  in  relation  to  those  of  uninvaded 
and  lightly  invaded  sites  (Fig.  4C),  but  PerMANOVA  found  invasion  level  to  be  only 
marginally  significant  when  soil  texture  and  habitat  were  included  in  the  model 
(Table  2). 


Within- site  evaluation  of  bacterial  community  variation  at  the  scale  of  plant 
neighborhoods  did  not  detect  any  significant  effects  of  plant  species  on  soil  bacterial 
communities  (Table  3).  For  fungi,  significant  PLANT-related  differences  in  community 
composition  were  found  in  three  of  the  ten  tested  sites  (Table  3;  Fig.  5).  Pairwise 
contrasts  revealed  that  fungal  communities  in  the  neighborhood  of  L.  cuneata  were 
different  from  those  of  bulk  soil  and  from  the  native  plant  species  at  the  lightly  invaded 
site  at  Ft.  Benning  (Fig.  5A).  Fungal  communities  in  the  neighborhood  of  L.  cuneata 
were  statistically  different  from  those  of  native  Chamaecrista  fasciculata  at  the  heavily 
invaded  site  of  Ft.  Stewart,  but  fungal  communities  associated  with  L.  cuneata  were  not 
distinguishable  from  bulk  soil  or  from  exotic  Medicago  lupulina  at  this  site  (Fig.  5B). 
Fungal  communities  associated  with  L.  cuneata  were  different  from  those  of  bulk  soil  at 
the  heavily  invaded  site  at  Ft.  Leonard  Wood,  while  the  fungal  communities  associated 
with  exotic  Coronilla  varia  were  not  distinguishable  from  either  (Fig.  5C). 


Table  3.  Significance  levels  of  small-scale  PLANT  effects3 


p-value 

INSTALLATION 

INVASION 

LEVEL 

Bacteria 

Fungi 

Stewart 

light 

heavy 

0.50 

0.17 

0.80 

0.02  * 

Benning 

light 

heavy 

0.17 

0.72 

0.008  ** 

0.24 

Campbell 

light 

heavy 

0.70 

0.56 

0.33 

0.15 

Leonard  Wood 

light 

heavy 

0.74 

0.83 

0.52 

0.04  * 

Riley 

light 

heavy 

0.48 

0.20 

0.40 

0.39 

“Significance  level  determined  by  permutational  MANOVA 
*  0.01  <  p  <  0.05;  **  0.001  <  p  <  0.01 


4.  Discussion 


4.1  Bacteria  respond  on  the  site-scale 

When  considering  sites  from  across  the  U.S.  invaded  range  of  L.  cuneata ,  we 
found  that  significant  changes  in  the  community  composition  of  soil  bacteria  were 
correlated  with  the  invasion  level  of  L.  cuneata  at  each  site  (Table  2;  Fig.  2).  At  sites 
where  L.  cuneata  dominated  the  plant  community,  the  soil  bacterial  communities  were 
consistently  different  from  those  of  uninvaded  sites  (Fig.  2),  and  these  effects  were  found 
across  the  range,  with  heavily  invaded  sites  from  all  installations  (used  to  block  the 
study)  being  shifted  in  the  same  direction  on  NMDS  plots  regardless  of  their  installation 
of  origin  (Figs.  2,  3).  Thus,  soil  bacteria  had  a  consistent  response  to  predominant  ground 
cover  of  L.  cuneata  in  spite  of  the  broad-scale  variation  in  the  underlying  microbial 
species  pool  that  may  result  from  other  unmeasured  sources  of  variation  (Table  1). 
Bacterial  communities  did  not  show  a  significant  response  to  either  the  SOIL  texture  of 
sites  nor  to  the  predominant  habitat  type  (Table  2),  and  the  response  to  L.  cuneata 
invasion  level  remained  significant  even  if  we  conservatively  account  for  these 
potential  sources  of  variability  prior  to  analysis  (Table  2).  The  consistency  of  these  shifts 
at  each  of  the  installations  sampled  strengthens  the  linkage  between  bacterial  community 
structure  and  L.  cuneata  dominance  of  the  plant  community,  and  this  highlights  the 
importance  of  studying  invader-microbe  dynamics  at  multiple  geographic  locations 
across  the  invader’s  range. 

The  influence  of  plant  species  dominance  in  determining  the  composition  of  the 
soil  community  has  both  theoretical  and  empirical  support  (Bever  1994,  Bever  2003, 
Mangan,  et  al.  2010a,  Mangan,  et  al.  2010b,  Kulmatiski,  et  al.  2011).  Bever  (2003) 
presents  a  coupled  plant  competition-soil  microbe  model,  where  each  plant  has  a 
particular  influence  on  microbial  community  composition,  and  the  equilibrium  microbial 
community  depends  on  the  relative  abundance  of  each  plant  in  the  community.  In  this 
model,  each  plant  can  interact  with  its  own  microbes  as  well  as  those  that  have  been 
selected  by  its  neighbors  (Bever  2003).  This  line  of  thinking  led  us  to  visualize  a  plant's 
influence  "bleeding  out"  from  its  root  zone  to  where  it  may  potentially  impact  the 
microbiota  and  other  plants  within  its  soil  "neighborhood."  While  plants  are  known  to 
interact  strongly  with  the  rhizosphere  communities  in  direct  contact  with  its  roots  (Kent 
and  Triplett  2002,  Berg  and  Smalla  2009),  it  is  on  the  neighborhood  scale  of  interaction 
that  changes  imposed  by  an  individual  plant  can  most  likely  influence  its  plant 
competitors.  We  examined  plant  neighborhood  effects  by  evaluating  microbial 
communities  from  the  root  balls  of  multiple  plant  species,  but  we  did  not  find  evidence 
that  bacterial  community  composition  varied  at  the  scale  of  plant  neighborhoods  (Table 
3).  Instead,  each  sample  derived  from  a  given  site  appeared  to  be  a  random  sample  of  the 
overall  bacterial  community  at  the  entire  site.  Therefore,  it  appears  unlikely  that  L. 
cuneata  invasions  are  supported  by  neighborhood-scale  bacterial  interactions  that 
influence  plant  competition.  However,  we  did  not  look  specifically  at  the  soil  in  direct 
contact  with  plant  roots,  i.e.  the  rhizosphere.  Thus,  we  cannot  say  whether  or  not  the  plant 
species  in  this  study  induced  bigger  changes  to  the  bacterial  communities  that  were  in 
direct  contact  or  near-direct  contact  with  root  surfaces.  Instead,  we  found  that  a 


significant  component  of  microbial  species  turnover  occurred  on  the  scale  of  entire  sites 
(-100-1000  m2). 

The  build  up  of  L.  cuneata  litter  at  heavily  invaded  sites  is  one  potential 
mechanism  by  which  this  plant  could  modify  soil  communities  on  the  scale  of  sites  but 
not  at  the  root  neighborhood  scale.  L.  cuneata  litter  contains  a  high  concentration  of 
phenolic  tannins  (Kalburtji,  et  al.  1999)  that  are  reported  to  have  bioactive  and 
antimicrobial  properties  (Langdale  and  Giddens  1967,  Kalburtji,  et  al.  2001,  Min,  et  al. 
2008).  The  build  up  of  these  compounds  in  soils  supporting  high  densities  of  L.  cuneata 
may  thus  directly  alter  the  soil  communities  by  changing  the  chemistry  of  soils,  as  has 
been  shown  for  several  other  invasive  plants  (Callaway,  et  al.  2004,  Callaway,  et  al. 
2008).  Kardol  and  colleagues  (2010)  report  that  stands  with  high  densities  of  L.  cuneata 
have  a  reduced  abundance  of  bacterivorous  soil  nematodes,  which  may  potentially  alter 
microbial  communities  via  food  web  effects.  Changes  in  the  relative  recalcitrance  of  litter 
can  also  affect  the  community  composition  of  soil  bacterial  communities  (Griffiths,  et  al. 
1999).  The  high  tannin  concentration  of  L.  cuneata  litter  has  been  shown  to  slow  down 
the  litter  decomposition  rate,  as  well  as  the  release  rates  of  N,  P,  K,  Mg,  and  other  trace 
elements  (Kalburtji,  et  al.  1999),  and  soils  in  L.  amrata-dominated  old  fields  have  been 
found  to  have  lower  soil  enzyme  activity  in  comparison  to  stands  dominated  by  the 
perennial  grass  Festuca  pratensis  (Kardol,  et  al.  2010).  L.  cuneata  associates  with 
nitrogen-fixing  bacteria,  and  its  dominance  at  a  site  may  also  influence  soil  nitrogen 
levels  or  change  the  overall  "quality"  of  litter,  both  of  which  can  influence  soil 
communities  and  play  important  roles  in  plant  invasion  and  competition  (Heneghan,  et  al. 
2006,  Miki,  et  al.  2010).  Nitrogen  fixation,  in  and  of  itself,  is  not  a  novel  trait  at  any  of 
these  sites,  as  all  of  the  plants  in  this  study  were  legumes,  but  dominance  of  L.  cuneata 
could  potentially  increase  the  overall  nitrogen  input  rates  at  a  site  due  to  its  high  nitrogen 
fixation  potential  (Lynd  and  Ansman  1993).  At  present,  we  do  not  have  the  data  to 
evaluate  whether  these  biogeochemical  changes  were  present  in  the  soils  under 
consideration  here,  but  our  finding  of  bacterial  community  shifts  in  relation  to  L.  cuneata 
site  dominance  suggests  that  future  work  on  this  invader  should  consider  changes  on  this 
whole-site  scale. 

The  high  final  stress  on  NMDS  plots  (Figs.  2-5)  indicates  that  higher  dimensional 
solutions  are  needed  to  accurately  portray  the  relationships  of  all  of  the  sample  points. 
However,  in  these  higher  dimensional  solutions,  the  configuration  of  the  treatment  group 
centroids  was  not  substantially  changed  from  what  is  portrayed  here,  and  so  we  have 
presented  two-  and  three-dimensional  solutions  here  for  ease  of  interpretation.  We  base 
our  conclusions  on  the  overall  and  pairwise  PerMANOVAs,  which  utilized  the  original 
sample  distance  matrices  and  thus  were  not  affected  by  NMDS  stress  or  dimensionality 
of  the  data.  Nevertheless,  the  high  final  stress  is  indicative  of  the  complex  structure  of  the 
community  data,  which  were  highly  variable  (note  the  large  standard  deviations  on  Figs. 
2-5),  with  a  weak  overall  correlation  structure.  Because  ARISA  fingerprinting  is  a  whole- 
community  approach,  our  dataset  included  taxa  that  were  correlated  with  L.  cuneata 
invasion  and  those  that  were  not.  The  impact  of  invasion  level  on  the  bacterial 
communities  was  significant  (Table  2),  but  the  inclusion  of  uncorrelated  taxa  in  the 
dataset  contributed  a  great  deal  of  noise  and  complexity,  making  the  overall  effect  subtle 


(Fig.  2).  Isolation  of  such  subtle  effects  represents  a  major  challenge  in  microbial 
community  ecology,  but  identification  of  the  proper  scale  of  response  can  be  a  useful 
guide  for  designing  more  powerful  experiments. 

An  important  caveat  of  our  work  is  that  the  analysis  of  this  observational  dataset 
does  not  identify  the  causal  direction  of  the  link  between  soil  bacterial  community 
composition  and  L.  cuneata  invasion  intensity.  Several  possibilities  can  account  for  the 
patterns  seen  here.  Soil  bacteria  may  be  responding  directly  to  changes  caused  by  L. 
cuneata  invasion,  perhaps  through  one  or  more  of  the  mechanisms  discussed  above. 
Alternatively,  L.  cuneata  may  reach  higher  densities  at  sites  with  a  particular  assemblage 
of  bacterial  taxa.  It  is  also  possible  that  both  L.  cuneata  and  the  soil  bacteria  may  be 
responding  to  some  unmeasured  source  of  site-level  variation  (e.g.  disturbance  history, 
soil  moisture),  but  the  consistency  of  the  bacterial  response  (Fig.  2)  across  such  a  large 
geographic  region  suggest  that  this  is  unlikely.  Further  experimental  work  will  be  needed 
to  evaluate  the  precise  mechanisms  at  play  and  to  determine  if  any  such  bacterial  changes 
on  the  whole-site  scale  have  lasting  consequences  for  L.  cuneata  fitness  or  the  fitness  of 
its  native  competitors. 

4.2  Fungi  respond  on  the  neighborhood  scale 

While  a  significant  component  of  bacterial  community  variation  could  be  linked 
to  L.  cuneata  invasion,  we  did  not  find  sufficient  evidence  to  link  soil  fungi  to  L.  cuneata 
invasion  level  at  the  scale  of  entire  sites  (Table  2).  This  is  primarily  because,  in 
contrast  to  bacteria,  fungal  communities  showed  significant  responses  to  soil  texture 
and/or  habitat  type  (Table  2;  Fig.  4),  both  of  which  were  collinear  with  invasion  level 
(Table  1).  To  handle  this  collinearity,  we  took  the  conservative  approach  of  partialling 
out  all  variation  related  to  soil  texture  and  habitat  type  before  conducting  the 
invasion  level  test  on  the  residuals.  In  contrast  to  bacteria,  invasion  level  only 
marginally  affected  fungal  communities  after  accounting  for  these  other  sources  of  site- 
level  variability  in  our  statistical  model.  The  soil  environment  and  the  predominant 
habitat  type  are  known  to  affect  mycorrhizal  fungal  communities  (Allen,  et  al.  1995, 
Smith  and  Read  2002).  In  our  study,  fungal  response  to  soil  texture  largely  reflected 
changes  along  a  sand-to-loam  gradient,  with  silty  loam  soils  (most  soils  in  this  study  had 
this  texture;  Table  1)  being  intermediate  to  these  groups.  Soil  fungal  communities  from 
forested  habitats  were  found  to  be  significantly  different  from  those  of  prairies  (Fig.  4A). 
Soil  texture  and  habitat,  in  conjunction  with  regional  variation  related  to  our 
blocked  design,  accounted  for  significant  variability  in  soil  fungal  communities,  and  they 
represent  an  alternative  hypothesis  to  invasion  level.  While  soil  fungi  may  be 
responding  to  L.  cuneata  densities  (marginally  significant  p-value,  Table  2),  our  evidence 
for  this  is  inconclusive.  Site-level  differences  in  fungal  community  composition  are 
equally  like  to  result  from  some  combination  of  geographic,  soil  textural,  or  habitat 
factors,  and  future  work  should  focus  on  teasing  apart  these  drivers. 

We  did  find  evidence  of  fungal  community  variation  at  the  plant  neighborhood 
scale  (Table  3;  Fig.  5).  Because  fungal  hyphae  are  larger  than  bacteria,  it  is  possible  that 
our  root  ball  sampling  included  some  fungi  extending  from  the  roots  of  the  plants,  and 


thus  fungal  rhizosphere  communities  are  better  represented  in  the  plant  soil  neighborhood 
than  are  bacterial  rhizosphere  communities.  If  soil  fungi  are  more  likely  than  bacteria  to 
respond  on  this  neighborhood  scale,  we  suggest  that  the  kinds  of  pairwise  plant-soil 
interactions  described  by  Bever  (2003)  may  be  more  likely  to  be  mediated  by  fungal 
networks  extending  from  one  plant  to  another  instead  of  by  localized  alterations  to  the 
soil  bacteria.  This  could  help  explain  why  significant  changes  to  fungal  communities 
have  been  more  consistently  identified  in  plant-feedback  studies  than  have  changes  in 
bacterial  communities  (Bezemer,  et  al.  2006a,  Bezemer,  et  al.  2006b,  Kardol,  et  al.  2007, 
van  der  Putten,  et  al.  2007b,  Lorenzo,  et  al.  2010).  However,  functionally-oriented 
feedback  studies  are  necessary  to  rigorously  test  this  hypothesis  for  L.  cuneata. 

We  found  significant  differences  between  the  fungal  communities  of  L.  cuneata 
and  those  of  bulk  soil  at  two  sites  (Fig.  5A,  C),  but  this  was  not  true  at  a  third  site  (Fig. 
5B).  Although  L.  cuneata  invasion  may  be  facilitated  by  beneficial  relationships  with 
mycorrhizal  fungi  (Lynd  and  Ansman  1993),  we  were  not  able  detect  these  host- 
associated  fungal  community  shifts  consistently  across  its  invaded  range.  We  did  find 
differences  between  the  neighborhoods  of  L.  cuneata  and  those  of  co-occurring  native 
plants  (Fig.  5A-B),  but  L.  cuneata  fungal  communities  were  not  deemed  to  be 
significantly  different  than  those  of  co-occurring  exotics  (Fig.  5B-C). 

Introduced  plants  represent  novel  foci  of  interaction  for  soil  microbial 
communities,  and  it  may  be  expected  that  novel  communities  would  develop  in  their 
vicinity.  Over  time,  a  history  of  repeated  interaction  between  plants  and  soil  fungi  can 
lead  to  shifts  in  the  fungal  community  composition  (Lankau  2011)  and  the  strength  of 
interaction  (Klironomos  2003).  Differences  in  "native"  vs.  "exotic"  fungal  communities 
may  play  an  important  role  in  the  invasion  process  if  exotic  plants  are  less  susceptible  to 
naive  fungal  pathogens  or  form  novel  mutualisms  with  mycorrhizal  fungi  (Richardson,  et 
al.  2000).  However,  without  functional  data,  we  are  not  able  to  determine  whether 
"native"  vs.  "exotic"  fungal  communities  have  important  implications  for  L.  cuneata 
invasion  dynamics. 

4.3  Conclusions 

The  implicit  consideration  of  multiple  spatial  scales  has  allowed  us  to  identify  the 
site  scale  (-100-1000  nT)  as  a  likely  scale  of  interaction  between  L.  cuneata  invasion  and 
soil  bacterial  communities,  while  fungal  communities  may  more  often  respond  at  the 
scale  of  localized  plant  "neighborhoods."  This  differential  scaling  of  bacterial  and  fungal 
association  with  L.  cuneata  invasion  may  indicate  that  soil  bacteria  and  soil  fungi  play 
different  ecological  roles  in  relation  to  this  invader.  Because  the  invasion  intensity 
relationship  with  soil  bacteria  was  evident  in  all  samples  from  heavily  invaded  sites 
regardless  of  plant  species,  we  propose  that  dominance  of  L.  cuneata  at  a  site  may  alter 
soil  bacterial  species  pools  at  sites  through  indirect  channels,  and  that  these  changes  may 
take  some  time  manifest.  On  the  other  hand,  fungi  were  more  likely  to  respond  at  the 
localized  scale  of  plant  neighborhoods,  and  this  scale  of  response  is  a  potential  indicator 
of  direct  plant-fungus  interactions  such  as  mutualism  and  parasitism.  Future  work  on  L. 


cuneata  invasions  should  take  into  consideration  these  different  scales  when  seeking 
potential  microbial  interactions  with  this  invader. 
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Figure  Legends 


Figure  1.  Typical  site  appearance  of  main  treatments.  These  images  show  typical 
plant  cover  and  habitat  conditions  under  different  invasion  levels  of  L.  cuneata.  (A) 
Uninvaded  forest  clearing  from  Ft.  Leonard  Wood,  Missouri;  (B)  lightly  invaded  open 
field  from  Ft.  Stewart,  Georgia  with  ellipses  indicating  individual  L.  cuneata  plants;  (C) 
heavily  invaded  open  field  from  Ft.  Benning,  Georgia  with  ellipses  indicating  several 
large  patches  of  L.  cuneata.  Almost  all  of  the  greenery  in  the  foreground  of  panel  (C)  is 
L.  cuneata.  In  contrast,  note  the  greater  plant  species  richness  and  evenness  apparent  in 
panels  (A)  and  (B).  Habitat  conditions  follow  the  definitions  used  in  Table  1. 

Figure  2.  Bacterial  community  composition  varies  with  L.  cuneata  invasion  level. 

Two-dimensional  nonmetric  multidimensional  scaling  (NMDS)  plots  reveal  how 
different  invasion  levels  of  L.  cuneata  relate  to  changes  in  soil  bacterial  communities. 
The  large  symbols  represent  the  centroids  of  all  samples  from  each  invasion  level 
group,  and  the  bars  show  the  standard  deviation  along  each  NMDS  axis.  Groups  with 
centroids  labeled  with  a  common  letter  were  not  deemed  to  be  different  in  pairwise  post- 
hoc  tests  (permutational  MANOVA)  using  the  Bonferroni  correction  to  maintain  an 
overall  a  =  0.05.  The  final  stress  is  reported  as  a  percentage. 

Figure  3.  Bacterial  and  fungal  communities  show  different  responses  to  geographic 
variation.  Three-dimensional  nonmetric  multidimensional  scaling  (NMDS)  plots  reveal 
geographic  structuring  in  soil  communities  of  (A,  B)  bacteria  and  (C,  D)  fungi.  The  large 
symbols  represent  the  centroids  of  all  samples  from  each  military  installation  (the 
Regional  variable),  and  the  bars  show  the  standard  deviation  along  each  NMDS  axis. 
Groups  with  centroids  labeled  with  a  common  letter  were  not  deemed  to  be  different  in 
pairwise  post-hoc  tests  (permutational  MANOVA)  using  the  Bonferroni  correction  to 
maintain  an  overall  a  =  0.05.  Note  that  panels  (A)  and  (B)  are  part  of  the  same  3- 
dimensional  configuration  and  share  a  common  stress  value  and  vertical  axis,  as  do 
panels  (C)  and  (D).  The  configurations  of  panels  (A)  and  (C)  are  qualitatively  similar  to 
those  resulting  from  two-dimensional  NMDS,  and  can  thus  be  compared  with  Fig.  2A 
and  2B,  respectively.  The  final  stress  is  reported  as  a  percentage. 

Figure  4.  Fungal  communities  vary  with  habitat  and  soil  type.  Two-dimensional 
nonmetric  multidimensional  scaling  (NMDS)  plots  reveal  how  soil  fungal  communities 
vary  with  respect  to  (A)  habitat,  (B)  soil  texture,  and  (C)  invasion  level  of  L. 
cuneata.  The  large  symbols  represent  the  centroids  of  all  samples  from  each  of  the 
respective  groups,  and  the  bars  show  the  standard  deviation  along  each  NMDS  axis. 
Groups  with  centroids  labeled  with  a  common  letter  were  not  deemed  to  be  different  in 
pairwise  post-hoc  tests  (permutational  MANOVA)  using  the  Bonferroni  correction  to 
maintain  an  overall  a  =  0.05.  Note  that  these  configurations  resulted  from  the  same 
NMDS  analysis,  and  the  stress  is  the  same  for  all  of  these  panels.  The  final  stress  is 
reported  as  a  percentage. 

Figure  5.  Some  plant  "neighborhoods"  influence  soil  fungal  communities.  Two- 
dimensional  nonmetric  multidimensional  scaling  (NMDS)  plots  reveal  small-scale 


differences  of  soil  fungal  communities  at  (A)  the  lightly  invaded  site  of  Ft.  Benning,  (B) 
the  heavily  invaded  site  of  Ft.  Stewart,  and  (C)  the  heavily  invaded  site  of  Ft.  Leonard 
Wood.  The  large  symbols  represent  the  centroids  of  all  samples  from  each  invasion 
level  group,  and  the  bars  show  the  standard  deviation  along  each  NMDS  axis.  Groups 
with  centroids  labeled  with  a  common  letter  were  not  deemed  to  be  different  in  pairwise 
post-hoc  tests  (permutational  MANOVA)  using  the  Bonferroni  correction  to  maintain  an 
overall  a  =  0.05.  The  final  stress  is  reported  as  a  percentage.  Abbreviations  of  plant 
species:  C.  var=  Coronilla  varia,  C.  fas  =  Charnaecrista  fasciculata,  L.  cun  =  Lespedeza 
cunecita,  L.  vir  =  L.  virginica,  M.  lup  =  Medicago  lupulina. 
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